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Abstract - An analysis ol the residual-velocity field of OB associ- 
ations within 3 kpc of the Sun has revealed periodic variations in the 
radial residual velocities along the Galactic radius vector with a typical 
scale length of A = 2.0 ± 0.2 kpc and a mean amplitude of /_r = 7 ± 1 
km s^^. The fact that the radial residual velocities of almost all OB- 
associations in rich stellar-gas complexes are directed toward the Galactic 
center suggests that the solar neighborhood under consideration is within 
the corotation radius. The azimuthal-velocity field exhibits a distinct 
periodic pattern in the < Z < 180° region, where the mean azimuthal- 
velocity amplitude is /g = 6 ± 2 km s~^. There is no periodic pattern of 
the azimuthal-velocity field in the 180 < I < 360° region. The locations 
of the Cygnus arm, as well as the Perseus arm, inferred from an analysis 
of the radial- and azimuthal-velocity fields coincide. The periodic pat- 
terns of the residual- velocity fields of Cepheids and OB associations share 
many common features. Key words: star clusters and associations, stellar 
dynamics, kinematics; Galaxy (Milky Way), spiral pattern. 

1. INTRODUCTION 

The location of the spiral arms in our Galaxy and their influence on the kinemat- 
ics of gas and young stars is undoubtedly of great importance for understanding 
the large-scale hydrodynamic processes and the evolution of stellar groupings 
and that of the Galaxy as a whole. Only in our own Galaxy it is possible to de- 
rive the space velocity field of young stars and analyze the radial and azimuthal 
velocity components simultaneously. However, even now the location of spiral 
arms in the Galaxy remains a subject of discussion. There are two main ap- 
proaches to the problem. The first one consists in identifying spiral arms from 
regions of enhanced density of young objects in the galactic disk (Morgan et 
al. 1952; Fenkart and BinggeU 1979; Humphreys 1979; Efremov 7198; Berd- 
nikov and Chernin 1999, and others). The main idea of the second approach is 
to look for telltale kinematical signatures of spiral arms (Burton 1971; Burton 
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and Bania 1974; Crczc and Mennessier 1973; Gcrasimcnko 1993; Mishurov et 
al. 1979, 1997; Mel'nik et al. 1998, 1999; Sitnik and Mel'nik 1999; Sitnik et 
al. 2001, and others). We consider the latter approach to be more promising, 
because kincmatical analyses can be performed on incomplete samples, whereas 
when comparing stellar space densities in different regions of the Galaxy, it is 
necessary to allow for observational selection, which is a very difficult problem. 

However, when performing kincmatical analyses, one is up against another 
problem. Almost all stars with known line-of-sight velocities and proper mo- 
tions are located within 3 kpc of the Sun, with the kincmatical data becoming 
extremely scarce at greater heliocentric distances. As a result, it is impossible 
to infer the pitch angle of spiral arms directly from observations, because such 
a determination would require kinematic data for a solar neighborhood com- 
parable in size with the distance to the Galactic center. The basic idea of our 
approach to analyzing the velocity field is to determine the wavelength A of pe- 
riodic velocity variations along the galactic radius-vector ignoring the nonzero 
spiral-arm pitch angle and assuming spiral-arm fragments to have the shapes of 
circular segments. The wavelength in question is to a first approximation equal 
to the interarm distance. We first applied this technique to the Cepheid velocity 
field and foimd A = 1.9 ±0.2 kpc (Mel'nik et al. 1999). It is this parameter and 
not the spiral-arm pitch angle that we determine directly from an analysis of 
stellar kinematics. To estimate the pitch angle, we must also know the number 
or arms m. We can then determine the mean pitch angle i of spiral arms in 
terms of a model of regular spiral pattern using the simple relation tan i = 
of the density wave theory (Lin et al. 1969). The number of spiral arms is very 
difficult to establish even from radio observations, because no reliable distance 
estimates are usually available for gas clouds and other spiral-arm indicators 
observed at radio frequencies at large heliocentric distances. For a two-armed 
spiral pattern a wavelength of A = 1.9 kpc corresponds to a mean pitch angle 
of i = 5°. 

In the density- wave theory there is a fixed phase shift of 7r/2 between the 

oscillations of radial and azimuthal components of velocity perturbations for 
tightly wound arms (Lin et al. 1969). However, simulations of the kinematics of 
spiral arms with allowance for gas-dynamic effects (Roberts 1969) showed both 
components to exhibit maximum variations at the shock front. Therefore, one 
would expect that the oscillations of the radial and azimuthal components of 
residual velocity must be synchronized, and their phases should not be shifted 
relative to each other. The allowance for shock effects, even in a coarse approx- 
imation, requires an independent determination of the phases of oscillations of 
the radial and azimuthal components of residual velocity. It is the dropping 
of the requirement of the fixed phase shift between the oscillations of the two 
components that allowed considerable perturbations in the azimuthal residual 
velocities to be found in the < I < 180° region. 

Unfortunately, our method allows us to determine the arm location only 
up to A/2, i.e., up to shifting arms into the interarm space. Moreover, it is 
impossible to choose between the two solutions for the location of the spiral 
pattern based on kincmatical data alone. For the final choice, we invoke the 
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additional information about the location of starburst regions with respect to 
the periodic pattern found in the velocity field studied. This information allows 
us not only to choose the right spiral-pattern solution, but also to determine 
the position of the region considered (within 3 kpc from the Sun) with respect 
to the corotation radius. 

It is impossible to analyze the spiral pattern of our Galaxy without the 
knowledge of reliable distances to the objects studied. The use of OB asso- 
ciations instead of individual stars allows, due to averaging, a more reliable 
distance scale to be constructed and more reliable velocities to be derived. OB 
associations arc sparse groupings of young stars (Ambartsumyan 1949). A com- 
parison of virial mass estimates of OB-associations with the masses estimated by 
modeling their stellar content suggests that these groupings are gravitationally 
unbound. Presently, there arc several partitions of galactic OB-stars into asso- 
ciations, those of Blaha and Humphreys (1989), Garmany and Stencel (1992), 
and Mel'nik and Efremov (1995). All these partitions are based on the cat- 
alog of luminous stars of Blaha and Humphreys (1980). However, Garmany 
and Stencel (1992) identified their OB-associations only in the 50 < ^ < 155° 
region. Mel'nik and Efremov (1995) used cluster analysis technique to identify 
the densest and most compact groups, the cores of OB-associations. However, 
these groupings contain twice as few stars than the associations of Blaha and 
Humphreys and, moreover, in dense regions, kinematical data are available for 
a smaller fraction of stars. We therefore consider the partition of OB-stars into 
associations suggested by Blaha and Humphreys (1989) to be more suitable for 
kinematical analyses. The sky-plane sizes of most of the OB-associations of 
Blaha and Humprcys (1989) do not exceed 300 pc (except Cep OBI and NGC 
2430), and the use of these objects for identifying periodic patterns with typical 
scale lengths greater than 1 kpc is a quite correct procedure. The inclusion of a 
list OB-stars in the HIPPARCOS (1997) program allowed their space motions 
to be analyzed for the first time (de Zeeuw et aL1999). 

2. OBSERVATIONAL DATA 

To construct the velocity field of OB associations, we used the following data: 

1. The catalog of stars in OB associations by Blaha and Humphreys (1989) 
and heliocentric distances rsH of these associations; 

2. The catalog of classical Cepheids (Berdnikov 1987; Bcrdnikov et al. 2000) 
with distances on the so-called short distance scale (Berdnikov and Efre- 
mov 1985); 

3. The solar Galactocentric distance i?o = 7.1 ± 0.5 kpc (Rastorguev et al. 
1994; Dambis et al. 1995; Glushkova et al. 1998); 

4. The short distance scale for OB associations, r = CSr^ij, which is consis- 
tent with the short distance scale for Cepheids (Sitnik and Mel'nik 1996; 
Dambis et al. 2001); 
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5. Stellar linc-of sight velocities from the catalog by Barbier-Brossat and 
Figon (2000); 

6. Proper motions of stars adopted from the HIPPARCOS (1997) catalog. 

The distances to OB associations that we use here and the solar Galactocen- 
tric distance were obtained on the same distance scale, which agrees with the 
short distance scale for Cepheids. 

The catalog of Blaha and Humphreys (1989) includes a number of extended 
open clusters, such as CoUinder 121 and Trumpler 16, because it is very diffi- 
cult to distinguish unambiguously between OB-associations and young clusters. 
However, we did not use in our analysis other young clusters that were not in- 
cluded in the catalog of Blaha and Humphreys (1989), because of a considerable 
overlap between the lists of OB-association and young-cluster stars with known 
kinematical parameters (Glushkova 2000, private communication). 

We determined median line-of-sight velocities for a total of 70 OB asso- 
ciations containing at least two stars with known line-of-sight velocities, and 
median tangential velocities, for 62 associations containing at least two stars 
with known proper motions 1^. The velocity of each OB association is based, on 
the average, on 12 line-of-sight velocities and 11 proper motions of individual 
stars. We excluded the distant association Ara OBIB (r — 2.8 kpc) from our 
sample because of its large Vz velocity component, which exceeds 20 km s~^. 



3. AN APPROACH TO THE SOLUTION 



In the case of tightly wound spiral arms, the velocity field must exhibit variations 
of the value and direction of residual stellar velocities (i.e., velocities corrected 
for the Solar apex motion and galactic rotation) on Galactocentric distance. 

We now write the expressions for the perturbation of radial Vu and azimuthal 
Vg components of residual velocities in the form of periodic functions of the 
logarithm of Galactocentric distance R: 

Vn = fRsmi^H^) + ^n), (1) 
A Ho 

Ve = fesm{^H^) + Vo), (2) 

where fu and fg are the amplitudes of variations of velocity components Vr and 
Vg. Parameter A (in kpc) characterizes the wavelength of the periodic velocity 
variations along the galactic radius-vector. The angels ipR and ipg determine 
the phases of oscillations at the solar Galactocentric distance. Assuming that 
galactic arms have the shape of logarithmic spirals, we adopted a logarithmic 
dependence of the wave phase on Galactocentric distance, which degenerates 
into a linear function RQln{R/Ro) Ki R — Rq if {R — Rq)/Ro is small. 

To demonstrate that the periodic pattern in the field of residual velocities is 
independent of the adopted model of circular rotation, we found the parameters 

^the catalog is available at jhttp:/ /Infml. sai.msu.ru/ ~anna/page3.html| 
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of the periodic pattern jointly with those of differential circular galactic rotation 
and the components of solar velocity. We inferred all these quantities from a 
joint solution of Bottlinger equations (Kulikovskii 1985) for line-of-sight veloci- 
ties and velocity components Vi (Vi = 4.738 [km s~^ kpc"-'^ (arcscc yr""'^)"'^] 
/i;r, where fii is the proper- motion component along the galactic longitude) with 
allowance for perturbations induced by the density wave: 



Vr = —{—Uq cos Z cos 6 + vq sin Z cos 5 + wq sin 6) + 
+Ro%{R - Ro)smlcoab + 
+0.5%iR - Rof sin Z cos 6 - 

-/Hsin(^^^ In(^) + ifiR) cos(/ + 6*) cos b + 

A il5 

+f0 sm( — - — ln( — ) + ipe) sm{l + 6) cos 6; 
A Hq 



(3) 



Vi = —{uq sinZ + vq cos6) + 

+n[{R- Ro){RoCosl -rcosb) + 
+0.5^0 (i? - Ro)'^{Ro cosl-r cos b) - flor cos b + 

+fR sin(^^ ^'^^R^^'^ + + 

+fesm{^lnA+^s)cos{l + 0); (4) 
A ito 

Here, 9 is the azimuthal Galactocentric angle between the directions toward the 
star and the Sun; is the angular velocity of galactic rotation at the solar 
Galactocentric distance; fif, and JIq are the first and second derivatives with 
respect to Galactocentric distance taken at a distance of Rq; Uq, Vq, and Wq 
are the solar velocity components relative to the centroid of OB-associations in 
the directions of X, Y, and Z axes, respectively. The X-axis is directed away 
from the Galactic center, the F-axis is in the direction of galactic rotation, 
and the Z-axis points toward the North Galactic Pole. Velocity components Uq 
and Vq include the solar-velocity perturbation due to the spiral density wave. 
(In galactic astronomy the X-axis is traditionally directed toward the Galactic 
center and one should therefore compare —Uq and not Uq with the standard 
solar apex). We adopted wq = 7 km s~^ for the solar velocity component along 
the Z-coordinate (Kulikovskii, 1985; Rastorguev et al. 1999). 

To linearize equations (3) and (4) with respect to the oscillation phases ipn 
and ife, we rewrite the formulas for perturbations of velocity components Vr 
and Vg as follows: 

Vr = Ar sin(— ^ ln( — )) + Br cos(— ^ ln( — )) (5) 

A riQ A riQ 

Ve = A, sin(^ In A) + Be cos(^ ln( A)) (g) 

A ±12 A xto 
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The parameters fa, fg, (fin and (fig can then be found from the relations: 



R — -^R + -^fli 



/(9 ~ + ^t'l 



(7) 



tan((pij) = Br/Ar; 



(8) 



We computed the weight factors p in the equations for Vr and Vi as follows: 



Here (To is the dispersion of residual velocities of OB associations with respect 
to the adopted model of motion (without allowance for the triaxial shape of the 
velocity distribution); evr and e^i are the standard errors of measured stellar 
line-of-sight velocities and proper motions, respectively. We determined the 
dispersion using iterations technique, which yielded uo = 6.6 km s~^. 

We then applied the least-squares technique to find a joint solution of the 
system of equations (3) and (4), which are linear with respect to the parameters 
Wq, Vq, Qq, f^o, f^o, Ar, Br, Ag, and Bg, with weight factors (9) and (10) 
and fixed A (see p. 499 in the book by Press et al. (1987)). We estimate the 
wavelength A by minimizing function X^Wi which is equal to the sum of squares 
of the normalized velocity residuals. 

4. RESULTS 

4-1- Parameters of the Rotation Curve and Periodic Pattern Inferred from the 
Entire Sample of OB Associations 

Figure 1 shows a-s a function of A, based on a joint solution of the system of 
equations (3) and (4) for line-of-sight and tangential velocities of OB associations 
located within 3 kpc from the Sun. takes its minimum value at A = 2.0 kpc. 
The resulting amplitudes of radial and azimuthal velocity perturbations are 
equal to fR = 6.6 ± 1.4 and fo = 1.8 ± 1.4 km s~^, respectively. Table 1 gives 
the inferred values of all determined parameters: Uq, vq, fl'^, Qq, SIq, A, Ar, Br, 
Ag, and Bg, as well as fR, fg, ipR and ipg computed using formulas (7) and (8). 
The table also gives the standard errors of the above parameters, the number 
N of equations used, and the rms residual ctq. 

We then performed numerical simulations in order to estimate the standard 
error of the resulting A. To this end, we fixed the actual galactic coordinates of 
OB-associations and simulated normally distributed random errors in the helio- 
centric distances of OB associations with standard deviations equal to 10% of 
the true distance, and then used formulas (3) and (4) to compute a theoretical 
velocity field with allowance for the perturbations due to the density wave (we 
took all parameter values from Table 1). We then added to the theoretical ve- 
locities simulated normally distributed random errors with a standard deviation 



Pvr = {al + elX'^' 



(9) 



Pvi = {(tI + (4738£^ir)2)-V2 



(10) 
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Table 

tern, 

the 


1. Parameters 
and solar-motion 
line-of-sight velocities 


of the circular rotation law, periodic pat- 
components inferred from an analysis of 
and proper motions of OB associations 


N 


"o 
km s"'^ 


vo 

km km 


rig r2o A 

kpc~^ km s"'^ kpc~* km s~^ kpc~^ kpc 


132 


-7.5 
±0.9 


11.2 
±1.3 


-5.0 1.5 30.2 2.0 
±0.2 ±0.2 =0.8 ±0.2 



Ar 


Br 


Ae 


Be 


fR 


fe 


•PR 


•fie 




km s^^ 


km s^^ 


kin s ^ 


km s^ '- 


km s^ '- 


km s ^ 


dcg 


dcg 


km s~^ 


5.2 


4.0 


1.5 


-1.0 


6.6 


1.8 


38° 


-33° 


6.6 


±1.4 


±1.3 


±1.4 


±1.4 


±1.4 


±1.4 


±12° 


±48° 





of a, which includes the contribution of observational errors {a = l/pvr and 
(7 = 1/pvi, see formulas (9) and (10)). Wc determined the wavelength A for 
each simulated velocity field and found the inferred A values to be unbiased and 
to have a standard error of 0.2 kpc. 

We also explored the possibility of periodic patterns emerging accidentally 
in the velocity field of OB associations due to chance deviations of individual 
velocities from the circular rotation law. To this end, we simulated random 
errors in the velocities and heliocentric distances of OB associations and deter- 
mined A, fft,, and f$ for each simulated field. Numerical simulations showed 
that 30% of all A values fall within the wavelength interval 1 < A < 3 kpc that 
is of interest for us. It is in this A interval that random fluctuations of the fleld 
of circular velocities can be attributed to density-wave effects. The mean ampli- 
tudes fa and fg are equal to 3 km s~^, i.e., about twice the standard errors of 
the corresponding parameters inferred for the actual sample of OB associations 
(Table 1). However, the probability of a periodic pattern with an amplitude of 
fft > 6.6 km s~^ and A in the 1 < A < 3 kpc interval emerging accidentally is 
extremely low P < 1%. Therefore the hypothesis about the periodic pattern 
with an amplitude equal to //{ = 6.6 km s~^ emerging as a result of chance fluc- 
tuations in the velocities and heliocentric distances of OB associations can be 
rejected at a confidence level of 1 — P > 99%. However, chance oscillations in the 
fleld of azimuthal velocities with amplitudes fe > 1.8 km s~^ and wavelengths 
A in the 1 < A < 3 kpc interval appear rather frequently, in 25% of the cases, 
and therefore the periodic pattern found in the field of azimuthal velocities with 
f$ = 1.8 km s~^ can well be interpreted in terms of random fluctuations. 

4-2. Location of Spiral Arms in the Galactic Plane 

A gravitational potential perturbation that propagates in a rotating disk at a 
supersonic speed produces a shock front, which affects the kinematics of gas and 
young stars born in this gas (Roberts 1969). The ages of OB associations do 
not exceed 5 x 10^ yr and, therefore, the motions of OB-associations must be 
determined mainly by the velocities of their parent molecular clouds (Sitnik et 
al. 2001). Inside the corotation radius the shock front must coincide with the 
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maximum radial velocity of streaming motions toward the Galactic center and 
maximum azimuthal velocity in the direction opposite that of galactic rotation. 
Velocities vary in value and reverse their direction as one recedes from the 
perturbation front. A sinusoidal law gives a first, coarse approximation to this 
pattern. According to the adopted model, in the interarm space, the radial 
and azimuthal components of streaming motions must be directed away from 
the Galactic center and along the galactic rotation, respectively (see Fig. 2 in 
Mel'nik et al. 1999). 

Figure 2 shows the distribution of OB associations with known space veloc- 
ities and the full vectors of residual velocities projected on the galactic plane. 
We determined residual velocities as the differences between heliocentric veloc- 
ities and circular rotation velocities, which also includes the motion of the Sun 
toward the apex, computed with parameters Uq, vq, f^o, ^o, ^a, adopted from 
Table 1. Also shown are circular arcs corresponding to the maximum mean 
radial velocity Vr toward the Galactic center as defined by formula (1) and fn 
and (fiu adopted from Table 1. Table 2 gives the following parameters for 59 
associations with known space velocities: radial (Vr) and azimuthal (Vg) com- 
ponents of residual velocities; components of residual velocities along the 
z-coordinate; Galactocentric R and heliocentric r distances, and galactic coor- 
dinates / and b. To characterize the reliability of velocities and distances listed 
in Table 2, we also give the numbers and ni of association stars with known 
line-of-sight velocities and proper motions, respectively, and also the number N 
of members of the OB associations used to determine distance r. 

Let us assume that the region studied is inside the corotation radius. It 
then follows, in view of the small value of the pitch angle, that the arcs shown 
in Fig. 2 should coincide with the shock front and must be located near the 
minima of gravitational potential minimum (Roberts 1969). Given a partition 
of yoimg galactic objects into stellar-gas complexes (Efrcmov and Sitnik 1988), 
one can identify the star-forming regions through which the arms drawn in Fig. 
2 pass. The arm located closer to the Galactic center passes in quadrant I near 
the OB associations of the Cygnus stellar-gas complex (Cyg 0B3, OBI, 0B8, 
and 0B9) and in quadrant IV, through the OB associations and young clusters 
of the stellar-gas complex in the constellations of Carina, Crux, and Centaurus 
(Car OBI, 0B2, Cru OBI, Cen OBI, Coll 228, Tr 16, Hogg 16, NGC 3766, and 
NGC 5606). Hereafter we refer to this arc as the Cygnus-Carina arm. Another 
arm, which is farther from the Galactic center, passes in quadrant II near the OB 
associations of the stellar-gas complexes located in the constellations of Perseus, 
Casiopeia, and Cepheus (Per OBI, NGC 457, Gas 0B8, 0B7, 0B6, 0B5, 0B4, 
0B2, OBI, and Cep OBI). In quadrant HI neither stellar-gas complexes nor 
even simply rich OB associations can be found to lie along the extension of this 
arc, which we refer to as the Perseus arm (Fig. 2). Note that Perseus-Cassiopeia 
and Carina-Centaurus stellar-gas complexes are the richest ones in the sense of 
the number of luminous stars their associations contain (see, e.g.. Table 2). 

If we assume that the solar neighborhood considered is located outside the 
corotation radius, the shock front and the minimum of potential should coin- 
cide with the maximum velocity of streaming motions directed away from the 
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Galactic center and maximum velocity of azimutlial streaming motions in the 
direction of galactic rotation. The arms should then be shifted by A/2 rela- 
tive to the lines drawn in Fig. 2, putting the rich stellar-gas complexes of the 
Cygnus-Carina and Perseus arms into the intcrarm space. Such a pattern would 
be inconsistent with modern concepts of star formation and results of observa- 
tions of other galaxies, which indicate that starburst regions concentrate toward 
spiral arms (Elmegrecn 1979: Efremov 1989). 

The fact that radial residual velocities of almost all OB associations in rich 
stellar-gas complexes are directed toward the Galactic center indicates that the 
region studied is located inside the corotation radius. 

4-3. Kinematically Distinct Star-Forming Regions 

However, things are not all that straightforward. An analysis of the data in 
Table 2 showed that about 30% of rich OB-associations (containing more than 
30 members with known photometric parameters, iV > 30), exhibit kinematic 
signatures characteristic of the interarm space. In particular, they have their 
radial residual velocity components Vr directed away from the Galactic center. 
This is not a surprise, because star formation can also proceed in the interarm 
space (Elmegreen and Wang 1987). In the solar neighborhood (Fig. 2) two 
regions can be identified where most of OB associations have radial velocities Vji 
directed away from the Galactic center. These are the associations of the Local 
system located in quadrants II and III (Vela 0B2, Mon OBI, Coll 121, Ori OBI, 
and Per 0B2) and the stellar-gas complex projected onto the constellations of 
Sagittarius, Scutim, and Serpens (Sgr OBI, 0B7, 0B4, Ser OBI, 0B2, Set 0B2, 
and 0B3) (Efremov and Sitnik 1988). These regions are located in the interarm 
space of the pattern shown in Fig. 2. It is the alternation of star-forming 
regions with positive and negative radial residual velocities Vr that determines 
the periodic pattern of the field of radial velocities of OB-associations. 

Within 3 kpc from the Sun a total of five star-forming regions can be iden- 
tified where almost all associations have the same direction of the radial com- 
ponent Vr of residual velocity. The contours of these regions are shown in Fig. 
2. Table 3 gives for each such region its mean Galactocentric distance R, mean 
residual velocities of associations Vr and V0; the interval of coordinates I and 
r, and the names of OB-associations with known space velocities it contains. 

Table 3 shows a well-defined alternation of the directions of the mean radial 
velocity Vr of OB associations as a function of increasing Galactocentric dis- 
tance R. The periodic pattern is especially conspicuous in Fig. 3a, which shows 
the variation of the radial component of residual velocity of OB assrciations 
along the Galactocentric distance. Radial velocities of OB associations in the 
Carina-Centaurus (R = 6.5), Cygnus (R = 6.9 kpc), and Perseus-Cassiopeia 
{R = 8.4) complexes are directed mainly toward the Galactic center, whereas 
those in the Sagittarius-Scutum complex (i? = 5.6 kpc) and in a part of the 
Local system {R = 7.4 kpc) are directed away from the Galactic center. The 
velocities of other OB associations located outside the above complexes are, on 
the average, smaller in magnitude, also in agreement with the periodic pattern 
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Table 3. Average residual velocities of OB-associations in the star-forming regions 



Region 


R, kpc 


Vr, ^ 
km s 


Ve, ^ 
km s 


I, deg 


r, kpc 


Associations 


Sagittarius 


5.6 


+11 ±3 


-i±i 


8-23 


1.3-1.9 


Sgr OBI, OB7, OB4, 
Ser OBI, OB2 
Set OB2, OB3; 


Carina 


6.5 


-6 ±2 


+5 ±3 


286-315 


1.5-2.1 


Car OBI, OB2, 
Cru OBI, Cen OBI, 
Coll 228, Tr 16, 
Hogg 16, 
NGC 3766, 5606; 


Cygnus 


6.9 


-7 ±3 


-11 ± 2 


73-78 


1.0-1.8 


Cyg OBI, OB3, OB8, 
OB9; 


Local System 


7.4 


+6 ±3 


+1±3 


160-360 


0.3-0.6 


Per OB2, Mon OBI, 
Ori OBI, Vela OB2, 
Coll 121, 140; 


Perseus 


8.4 


-7 ±2 


-5 ±2 


104-135 


1.8-2.8 


Per OBI, NGC 457, 
Gas OB8, OB7, OB6, 
OB5, 0B4, OB2, OBI, 
Cep OBI; 



inferred. 

4-4- Specific Features of the Periodic Pattern in the Velocity Field of OB Asso- 
ciations in the I < 180° and I > 180° Regions 

To study the specific features of the velocity field of OB associations, we ana- 
lyzed residual velocities of these objects as a function of Galactocentric distance 
separately for the two regions / < 180° and / > 180° (Fig. 4 and 5, respectively). 
Figure 4b shows the azimuthal velocity field of OB associations to exhibit a well- 
defined periodic pattern in the region I < 180°, whereas no such pattern can 
be seen in the field of velocity components Vg of the entire sample of OB as- 
sociations (Fig. 3b and Table 1). The mean amplitude of azimuthal velocity 
variations in the region considered is as high as fg = 5.1 ± 1.7 km s~^, i.e., 
almost triple the value of fg = 1.8 ± 1.4 km s~^ inferred from the entire sam- 
ple of OB associations. One can see two well-defined minima at Galactocentric 
distances R = 7.0 and R = 8.4 kpc. It is evident from a comparison of Figs. 
4a and 4b that the minima in the distributions of radial and azimuthal residual 
velocities are located at the same Galactocentric distances. 

In the density-wave theory including shock effects, the minima in the dis- 
tributions of the radial and azimuthal components of residual velocities must 
coincide with the shock front and should be located in the vicinity of the line 
of minimum potential (Roberts 1969). The striking agreement between the 
positions of minima as inferred from the distributions of radial and azimuthal 
residual velocities of OB associations of quadrants I and II (Figs. 4a, 4b) can 
be explained by the a shock front. The positions of these minima determine the 
kinematical positions of the Cygnus and Perseus-Cassiopeia arms putting them 
at Galactocentric distances of i? = 6.8 — 7.0 and i? = 8.2 — 8.5 kpc, respectively. 
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Table 4. 


Parameters 


of the 


periodic 


pattern 


in the ve- 


locity field 


of 


OB-associations 


located 


in 


different regions 


Region 


AT 


X, kpc 


fR 


fe 














km s ^ 


km s ^ 


deg 


deg 


km 


0<l< 180° 


73 


1.7 


6.7 


5.1 


52 


-15 


6.2 






±0.2 


±1.7 


±1.7 


±15 


±20 




30<l< 180° 


56 


1.4 


6.9 


6.1 


-7 


-26 


6.0 






±0.2 


±1.8 


±1.9 


±15 


±18 




180 < i < 360° 


59 


2.4 


8.1 


4.4 


7 


251 


6.2 






±0.4 


±2.0 


±1.8 


±12 


±24 





In the region / > 180° the periodic pattern is represented by a single mini- 
mum and two maxima in the distribution of radial velocities Vr (Fig. 5a). The 
minimum at the Galactocentric distance R = 6.2 — 6.5 kpc determines the kine- 
matic position of the Carina arm, whereas another minimum is absent, which 
would correspond to the extension of the Perseus arm toward quadrant III. The 
maximum at i? = 7.5 kpc (Fig. 5a) is associated with the positive radial veloc- 
ities of the Local system OB associations in quagrant III. No periodic pattern 
in the field of azimuthal velocities can be seen in the region / > 180° (Fig. 5b). 
That is why merging the association samples from the two regions {I < 180° 
and I > 180°) washes out the periodic pattern (Fig. 3b), although the latter is 
clearly outlined by the associations of quadrants I and II (Fig. 4b). 

For a quantitative analysis, we inferred the parameters of the periodic pat- 
tern in the velocity field of OB associations in two regions: < Z < 180° and 
180 < I < 360°, by solving the system of equations (3) and (4) with weight 
factors (9) and (10) and the parameters of circular rotation and solar-motion 
components adopted from Table 1. Table 4 gives the following parameters of 
the periodic pattern inferred for the two regions: A, fn, fe, fR and ipg, their 
standard errors, the number N of equations used, and the mean residual cfq. 

It is evident from Table 4 that in the region < Z < 180° radial and 
azimuthal residual velocities have similar variation amplitudes equal to Jr = 
6.7± 1.7 and fg = 5.1 ±1.7 km s""*^, respectively. The phases of the variations of 
the radial and azimuthal velocity components {ipR = 52±15 and (p$ = — 15±20) 
differ significantly from each other, although, as is evident from Figs. 4ab, the 
minima of the radial and azimuthal velocities are located at the same Galac- 
tocentric distances. The discrepancy is primarily due to the simple sinusoidal 
law adopted for our analysis of the periodic pattern. Large radial velocities of 
the OB-associations in the Sagittarius-Scutum region {0 < I < 30°) break the 
almost ideal periodic pattern outlined by the objects located in the Cygnus and 
Perseus arms and in the adjoining interarm space. Excluding from our sample 
the OB-associations located in the region < I < 30° changes phase ipR signif- 
icantly. The parameters of the periodic pattern inferred in the 30 < Z < 180° 
sector for the objects located in the Cygnus and Perseus-Casseopeia arms and 
in the adjoining interarm space are also listed in Table 4. It is evident from 
this table that in the sector considered the phases of radial and azimuthal ve- 
locity oscillations agree with each other within the errors {(pR — — 7 ± 15 and 
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tfo = — 26± 18, respectively). The reliably determined wavelength A = 1.4 ±0.2 
kpc for the region 30 < I < 180° is equal to the distance between the Cygnus 
and Perseus arms, or rather to that between the minima in the distributions of 
both radial and azimuthal residual velocities (Figs. 4a. 4b). 

In the region 180 < I < 360° the mean amplitude of radial velocity oscilla- 
tions is equal to fn = 8.1 ± 2.0 km s~^. Here the wavelength is determined as 
the distance between the maxima in the distribution of radial velocities, i.e., A 
proves to be equal to the distance between the interarm-space objects (Fig. 5a). 
The wavelength inferred for this region, A = 2.4 ± 0.4 kpc, is rather uncertain, 
because the Galactocentric dependence of residual velocity differs appreciably 
from the sinusoidal law. 

Numerical simulations of the velocity field allow the hypothesis that the 
variations of either radial or azimuthal velocities in the region 30 < I < 180°, 
with amplitudes equal to /r = 6.9 and /e = 6.1 km s~^, respectively, are due to 
accidental errors in the velocities and distances, to be rejected at a confidence 
level of 1 — P > 95%. On the other hand, the only statistically significant 
periodic pattern in the region 180 < I < 360° is that of radial velocities {l — P> 
99%), whereas azimuthal velocity variations can well {P = 15%) be interpreted 
in terms of random fluctuations. 

Figure 6a illustrates the specific features of the periodic pattern found in 
this work. It also shows the field of residual velocities of OB-associations and 
the circular arcs corresponding to the minima of residual radial (Vr) (solid line) 
and azimuthal (Vg) (dashed line) velocities based on parameters A, fu, fe, ^Pr 
and ipe, for regions 180 < Z < 360° and 30<l< 180° (Table 4). The arcs must 
be located in the vicinity of the lines of minimum gravitational potential. In the 
region 30 < I < 180° these lines determine the loci of the Cygnus and Perseus- 
Cassiopeia arms and in the region 180 < I < 360°, that of the Carina-Centaurus 
arm. Numerical simulations showed the inferred Galactocentric distances of 
arms and, correspondingly, the radii of arcs in Fig. 6a, to have standard errors 
of 0.1 - 0.2 kpc. 

Figure 6a illustrates all three specific features of the velocity field of OB- 
associations. First, the periodic pattern of azimuthal velocities in the region 
< Z < 180° and the absence of such pattern in the region 180 < / < 360°. 
Second, the agreement of Galactocentric distances of the Cygnus and Perseus- 
Cassiopeia arms as inferred from analyses of the fields of radial and azimuthal 
velocities in the region 30 < Z < 180°. Third, a 0.3 kpc shift of the kinematical 
positions of the Carina arm [R = 6.5 ± 0.1 kpc) relative to that of the Cygnus 
arm [R = 6.8 ± 0.1 kpc). Whether this shift is statistically significant remains 
an open question. 

5. COMPARISON OF PERIODIC PATTERNS IN THE VELOCITY 
FIELDS OF CEPHEIDS AND OB ASSOCIATIONS 

An analysis of the velocity field of Cepheids (Mel'nik et al. 1999) revealed a 
periodic pattern along the galactic radius- vector with a scale length of A = 1.9 ± 
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Table 5. 


Parameters 


of 


the periodic 


pattern 


in the 


velocity field 


of 


Cepheids 


located 


in 


different 


regions 


Region 


N 


A, kpc 


Jr 


fe 






o-() 








km s ^ 


km s ^ 


deg 


deg 


km s"'^ 


0<l< 180° 


217 


1.8 


6.7 


5.1 


85 


13 


10.2 






±0.2 


±1.6 


±1.5 


±14 


±17 




30<l< 180° 


165 


1.8 


6.5 


6.5 


82 


9 


9.8 






±0.2 


±1.8 


±1.7 


±16 


±15 




180 < i < 360° 


208 


1.8 


5.8 


2.0 


78 


154 


10.4 






±0.3 


±1.7 


±1.5 


±16 


±48 





0.2 kpc and mean oscillation amplitudes of fu = 6.2 ± 1.2 and fg — 2.1± 1.2 km 
s~^. It would be interesting to sec whether the periodic pattern of the Cepheid 
velocity field exhibits the same specific features as that of OB-associations. 

To answer this question, we inferred the parameters of the periodic pattern 
of the Cepheid velocity field in three regions < / < 180°, 30 < I < 180°, and 
180 < I < 360° in the same way as we did it for OB-associations (Table 5). 
Note that the exclusion of the region < I < 30° has no effect on the inferred 
parameters of the periodic pattern in the Cepheid velocity field. We used the 
parameters of circular motion and solar-motion components that we inferred 
from the analysis of the entire sample of Cepheid located within 3 kpc from the 
Sun [see Table in Mel'nik et al. (1999)]: = -7.6±0.8 km s'^; vq = 11.6±1.0 
km s-i; fl'o = -5.1 ± 0.2 km s'^ kpc'^; fl'^ = 1.0 ± 0.2 km s"^ kpc'^, and 
O = 29 ± 1 km s~^ kpc~^. These parameters agree within the errors with the 
corresponding parameters of the velocity field of OB-associations (Table 1). 

It is evident from Table 5 that in the region 30 < ? < 180° both radial 
and azimuthal velocity fields contain a periodic component. The amplitudes 
of velocity variations are equal to Jr = fe = 6.5 ±1.8 km s~^. Numerical 
simulations show that the hypothesis about the accidental nature of the periodic 
variations found in the fields of radial and azimuthal residual velocities in the 
region 30 < Z < 180° can be rejected at a confidence level of 1 — P > 99%. The 
periodic pattern revealed in the radial-velocity field in the region 180 < I < 360° 
with an amplitude of Jr = 5.8 ± 1.7 km s~^ cannot be also due to random 
fluctuations (1 — P > 99%), however the variations of azimuthal velocities with 
an amplitude of fe = 2.0 ± 1.5 km s~^ can well (P = 25%) be due to random 
fluctuations. 

Figure 6b shows the field of Cepheid residual velocities and, based on A, /h, 
fe, ipR, and ipg from Table 5, the circular arcs corresponding to the minima of 
mean radial (Vr) (solid line) and azimuthal (Vg) (dashed line) velocities. These 
arcs determine the kinematical positions of the Cygnus, Perseus-Cassiopeia, and 
Carina-Centaurus arms. The radii of arcs in Fig. 6b are determined with an 
accuracy of 0.1 — 0.2 kpc. 

The fields of residual velocities of OB-associations and Cepheids can be 
seen to have many features in common. First, the galacticentric distances of 
the Carina-Centaurus {R = 6.3 — 6.5 kpc) and Perseus-Cassiopeia {R = 8.1 — 
8.2 kpc) arms as inferred from analyses of radial velocities of Cepheids and 
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OB-associations agree well with each other. Second, both Cephcids and OB- 
associations located in the quadrants I and II exhibit periodic variations of 
the azimuthal velocity with an amplitude of » 6 ± 2 km s~^, whereas the 
azimuthal velocities of both populations show no periodic pattern in quadrants 
III and IV. However, the velocity field of Cepheids somewhat differs from that 
OB associations. Thus Cepheids do not show the 0.3 kpc shift between the 
positions of the Cygnus and Carina arms as inferred from radial velocities of 
OB associations. The Galactocentric distances of the Cygnus arm as inferred 
from analyses of radial velocities of Cepheids (R = 6.3 ±0.1 kpc) and OB- 
associations (R = 6.8 ± 0.1 kpc) differ significantly from each other. In contrast 
to OB associations, Cepheids do not show such a good agreement between the 
positions of the Cygnus and Perseus-Cassiopeia arms as inferred from analyses 
of the radial and azimuthal velocity fields. 

CONCLUSION 

An analysis of the field of space velocities of OB associations located within 3 
kpc from the Sun revealed periodic variations in the magnitude and direction 
of radial residual velocities Vr along the galactic radius-vector with a typical 
scale length of A = 2.0 ± 0.2 kpc and a mean amplitude of Jr = 7 ± 1 km s~^. 

We revealed five kinematically distinct star-forming regions where almost all 
OB-associations have the same direction of radial residual velocity Vr. The ra- 
dial velocities of OB associations in the Carina-Centaurus, Cygnus, and Perseus- 
Cassiopeia regions are directed mainly toward the Galactic center, whereas those 
of the Sagittarius-Scutum complex and of a part of the Local system, are di- 
rected away from the Galactic center. It is the alternation of star-forming regions 
with positive and negative radial velocities Vr that determines the periodic pat- 
tern of the radial velocity field of OB associations. 

The fact that rich Carina-Centaurus and Perseus-Casiopeia stellar-gas com- 
plexes lie in the vicinity of the minima in the distribution of radial velocities of 
OB-associations indicates that the region considered is located inside the coro- 
tation radius. The enhanced density of high-luminosity stars in these regions 
cannot be due to observational selection, because the Carina-Centaurus and 
Perseus-Cassiopeia complexes arc the most distant ones among those considered 
in this paper (Fig. 2 and Table 3). Furthermore, these regions cover sky areas 
extending for several tens of degrees and at a heliocentric distance of 2 kpc the 
mean extinction averaged over such large sectors depends little on the direction 
within the galactic plane. There is no doubt, the enhanced density of high- 
luminosity stars in the Carina-Centaurus and Perseus-Cassiopeia complexes is 
real and not due to the extremely low extinction along the corresponding lines 
of sight. 

The fact that the Perseus-Cassiopeia complex {R = 8.4 kpc) is located inside 

the corotation radius imposes a lower limit on the corotation radius, | Rc~Ro |> 
1.3 kpc and an upper limit on the spiral pattern speed, flp < 25 km s~^ kpc~^. 
Our conclusion about the corotation radius Rc being located in the outer 
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part of the Galaxy beyond the Perseus-Cassiopeia arm is in conflict with the 
conclusions of Mishurov and Zenina (1999a, 1999b) who argue that the Sun 
is located near the corotation radius | Rq — Rc \< 1 kpc. The latter authors 
based their conclusions on the small value of radial ( ffj = 2 ± 1 km s"-'^) and 
large value of azimuthal {fe=8zLl km s~^) velocity perturbation amplitudes 
they inferred (Mishurov and Zenina 1999b). Our analysis yielded a reversed 
amplitude proportion with /^j = 7 ± 1 and = 2 ± 1 km s"'^ both for Cepheid 
(Mel'nik et al. 1999) and OB-association data (Table 1). At the same time, 
the pitch angle found by Mishurov and Zenina (1999a, 1999b) for a two-armed 
spiral pattern corresponds to A = 2 kpc and is consistent with our results. 

As for the periodic pattern in the field of azimuthal velocities of OB-associations, 
it is observed in the quadrants I and II and is absent in the quadrants III and 
IV. The mean amplitude of azimuthal velocity variations in the region of the 
Cygnus and Perseus arms (30 < I < 180°) is as high as fe = 6 ± 2 km s~^, 
i.e., triple that of the entire sample of OB associations. Another specific fea- 
ture is the striking agreement between the positions of the Cygnus and Perseus 
arms as inferred from separate analyses of radial and azimuthal velocity fields 
of OB associations. This feature can be explained by shocks that develop when 
a density wave propagates through gas at a supersonic velocity. 

The periodic patterns in the residual velocity fields of Cepheids and OB- 
associations have very much in common: similar scale length of radial velocity 
variations along the galactic radius-vector, A = 2 ± 0.2 kpc and similar am- 
plitudes of velocity variations. Moreover, in both cases the azimuthal velocity 
fields of objects located in quadrants I and II exhibit a periodic pattern, whereas 
no such pattern can be seen in the azimuthal velocity field in the quadrants III 
and IV. 

The kinematical positions of the Carina {R = 6.3 — 6.5 kpc) and Perseus 
{R = 8.1 — 8.2 kpc) arms can be confidently inferred from an analysis of the 
field of radial velocities of both OB associations and Cepheids. It is the distance 
between these two arm fragments that determines the scale length of the vari- 
ations of the radial velocity component along the galactic radius-vector, A = 2 
kpc. 

The wavelength value that we inferred, A = 2 kpc, seems to coincide with 

that of the most unstable mode of galactic disk oscillations at the given Galacto- 
centric distance. The very existence of the spiral pattern suggests that the galac- 
tic disk is marginally unstable at the solar Galactocentric distance. Presently, 
the gaseous component of the galactic disk is considered to be instrumental 
in maintaining such a marginal instability (Jog and Solomon 1984; Bertin and 
Romeo 1988; Bertin et al. 1989). Jog and Solomon (1984) showed that the wave- 
length of the most unstable disk mode depends on the gas fraction. They found 
that the wavelengths of the most unstable disk modes at the solar Galactocen- 
tric distance should lie in the 1 — 5 kpc interval, where 1 and 5 kpc correspond 
to a purely gaseous and purely stellar disk, respectively. Our result A = 2 
kpc suggests that both components play important part in the dynamics of our 
Galaxy. 
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Table 2. Residual velocities of OB associations 



Association / b r R Vn Vg Vz, rir ni N 

kpc kpc km km km 



Sgr OB5 


o.°o 


-l.°2 


2.4 


4.7 


7.5 


9.0 


8.0 


2 


3 


31 


Sgr OBI 


7.6 


-0.8 


1.3 


5.8 


7.3 


-5.4 


1.0 


37 


29 


66 


Sgr OB7 


10.7 


-1.6 


1.4 


5.7 


9.2 


3.9 


-14.4 


3 


2 


4 


Sgr OB4 


12.1 


-1.0 


1.9 


5.2 


4.8 


-0.0 


-0.1 


9 


3 


15 


Ser OBI 


16.7 


0.1 


1.5 


5.7 


13.0 


-0.9 


1.4 


17 


12 


43 


Set OB3 


17.3 


-0.7 


1.3 


5.8 


2.7 


1.7 


3.0 


8 


3 


10 


Ser OB2 


18.2 


1.6 


1.6 


5.6 


13.7 


-1.9 


-3.7 


7 


5 


18 


Set OB2 


23.2 


-0.5 


1.6 


5.7 


26.3 


-2.3 


4.9 


6 


6 


13 


Vul OBI 


60.4 


0.0 


1.6 


6.5 


6.4 


-3.3 


6.1 


9 


8 


28 


Vul OB4 


60.6 


-1.2 


0.8 


6.7 


2.7 


-0.1 


-0.1 


3 


3 


9 


Cyg 0B3 


72.8 


2.0 


1.8 


6.8 


-14.5 


-7.8 


-4.0 


30 


18 


42 


Cyg OBI 


75.8 


1.1 


1.5 


6.9 


-3.9 


-8.3 


1.2 


34 


14 


71 


Cyg OB9 


77.8 


1.8 


1.0 


7.0 


-5.8 


-12.2 


0.2 


10 


8 


32 


Cyg OB8 


77.9 


3.4 


1.8 


7.0 


-3.9 


-13.9 


12.8 


9 


10 


21 


Cyg OB4 


82.7 


-7.5 


0.8 


7.0 


15.5 


4.8 


2.7 


2 


2 


2 


Cyg OB7 


89.0 


0.0 


0.6 


7.1 


2.5 


2.4 


3.7 


21 


28 


29 


Lac OBI 


96.7 


-17.7 


0.5 


7.2 


-3.6 


-2.8 


1.6 


2 


2 


2 


Cep OB2 


102.1 


4.6 


0.7 


7.3 


-3.2 


-0.1 


2.8 


37 


47 


59 


Cep OBI 


104.2 


-1.0 


2.8 


8.2 


-8.7 


-14.3 


3.1 


17 


24 


58 


Gas OB2 


112.0 


0.0 


2.1 


8.1 


-19.0 


-3.3 


6.2 


7 


5 


41 


Cep OB3 


110.4 


2.6 


0.7 


7.4 


-3.9 


-4.9 


0.5 


18 


15 


26 


Cas OB5 


116.1 


-0.5 


2.0 


8.2 


-12.7 


-2.7 


-11.2 


16 


13 


52 


Cas OB4 


120.1 


-0.3 


2.3 


8.5 


3.4 


0.8 


-6.4 


7 


7 


27 


Cas OB14 


120.4 


0.7 


0.9 


7.6 


11.6 


-2.7 


2.3 


4 


3 


8 


Cas OB7 


123.0 


1.2 


2.0 


8.4 


-10.3 


-9.2 


-2.7 


4 


8 


39 


Cas OBI 


124.7 


-1.7 


2.0 


8.4 


-6.8 


-2.5 


-7.0 


5 


3 


11 


NGC 457 


126.7 


-4.4 


2.0 


8.4 


-1.7 


1.0 


-6.7 


4 


2 


4 


Cas OB8 


129.2 


-1.1 


2.3 


8.7 


2.3 


2.5 


-3.7 


14 


9 


43 


Per OBI 


134.7 


-3.2 


1.8 


8.5 


-7.8 


-12.2 


-5.5 


81 


63 


167 


Cas OB6 


135.0 


0.8 


1.8 


8.4 


-11.2 


-6.8 


-3.8 


12 


13 


46 
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Fi gurc 1: The function X^iX) for the joint solution of the system of equations (3) and (4) 
for the entire sample of OB associations 
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Figure 2: The field of space velocities of OB-associations projected onto the galactic plane. 
The X-axis is directed away from the Galactic center, and the Y-axis is in the direction of 
galactic rotation. The Sun is at the origin. The circular arcs correspond to the maximum 
radial component Vfi of residual velocity toward the galactic center. One can sec five star- 
forming regions where almost all OB-associations have the same direction of radial velocity 
Vr. 
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Figure 3: 4, 5. Residual velocities Vn and Vg of OB-associatioiis as a function of Galacto- 
centric distance R. The whole sample and the regions < I < 180° and 180 < / < 360° are 
concidered. 
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Figure 6: The field of residual velocities of (a) OB-associations and (b) Cepheids. The 
circular arcs correspond to the positions of the arm fragments as inferred from analyses of 
radial (solid line) and azimuthal (dashed line) residual velocities in the regions 30 < / < 180° 
and 180 < 360°. 
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